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P^ ' We investigate chaotic behavior in a 2-D Hamiltonian system - oscihators with anharmonic 
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couphng. We compare the classical system with quantum system. Via the quantum action, we 
construct Poincare sections and compute Lyapunov exponents for the quantum system. We find 
that the quantum system is globally less chaotic than the classical system. We also observe with 
increasing energy the distribution of Lyapunov exponts approaching a Gaussian with a strong 
correlation between its mean value and energy. 
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I. INTRODUCTION 

Quantum chaos, which denotes properties of quantum systems whose underlying clas- 
sical system is chaotic, has been experimentally observed in irregular energy spectra of 
nuclei, of atoms perturbed by strong electromagnetic fields yy]. In a weak exterior mag- 
netic field, however, a hydrogen atom shows a level distribution being neither Poissonian 
nor Wignerian. How can we then compare classical chaos with quantum chaos? Is the 
quantum system more or less chaotic than the corresponding classical system? An under- 
standing of how classically regular and chaotic phase space is reflected in quantum systems 
is an open problem, since semiclassical methods of quantisation (EKB, Gutzwiller's trace 
formula) are not amenable to mixed dynamical systems ^. In order to describe nearest 
neighbor energy level spacing duistributions for mixed systems, workers have used different 
approaches. Percival p| has suggested a distribution based on a division of phase space 







into regular and irregular regions. Brody [^ has suggested a phenomenological fit between 

n 

a Wigner and a Poisson distribution. Izrailev |.5] suggested a distribution based on the con- 
cept of quantum localisation. Another interpolation between ensembles has been proposed 
by Lenz and Haake [6]. The interpolating distribution proposed by Berry and Robnik |2| 
uses an interpolation parameter determined from classical phase space. 

In order to address the above questions, workers have suggested concepts to reintroduce 
trajectories into quantum mechanics. One possibility is the use of the standard effective 
action, which takes into account quantum corrections of the classical action. Cametti et al. 
P] have considered the model of the 2-D anharmonic oscillator and computed the effective 
action via loop (h) expension to low order. Another approach is Bohm's interpretation 
of quantum mechanics which expresses the Schrodinger equation in polar form. It allows 
to introduce the concept of a particle trajectory and hence a quantum equivalent phase 
space. This approach has been used by Schwengelbeck and Faisal [9| to describe chaos in 
the kicked rotor and later by other workers for the 2-D anisotropic harmonic oscillators 
[ll|, the hydrogen atom in an external electromagnetic field [ll| and coupled anharmonic 



oscillators 12]. For the kicked rotor in 1-D, which is classically chaotic, it turns out that 
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uantum system is nonchaotic, i.e. the quantum Kolmogorov- Sinai entropy goes to 
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f rediscovered in 



12 1). In contrast to that, the 2-D anisotropic harmonic oscillator, 
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being classically non-chaotic, has been reported to exhibit quantum chaos [10||. 

A comparison with respect to chaotic behavior of the classical system and the quantum 
system has been investigated also in field theoretic models. Casetti et al. 1.3] considered 
the large N limit of A^-component 0^ oscillators in the presence of an external field. Using 
mean field theory they observed a strong suppression of quantum chaos due to quantum 
corrections moving the system away from a hyperbolic fixed point responsable for classical 
chaos. Matinyan and Miiller |lj] considered the model of massless scalar electrodynamics, 
which is classically chaotic. Using effective field theory and loop expansion, they noticed 
that quantum corrections increase the threshold for chaos by modifying the ground state 
of the system. 

To address the above questions, we apply the concept of the quantum action, de- 
veloped in Refs. 15|, |l6|, |l7|, |l8|, |l9l |20[. It expresses quantum transition amplitudes 



G{xfi, tfi] Xin, tin) in terms of a local action - the quantum action - S[x] defined by 
G{xfi, t = T] Xin, t = 0) = ^ Zexp[iS[xtraj]/h] , 

trajectories 

S[x]= / dt -mx^ - V{x{t)) , 
Jo 2 

Xtraj ■ SS[Xtraj] = . (l) 

Here 5* is evaluated along its trajectory (stationary point of S) going from boundary point 
{xin,t = 0) to {xfi,t = T). There r nay be several such stationary points. The quantum 
action has been explored numerically J15l.ll6l.ll7l.l20l| for confinement-type potentials (bound 



state spectrum), e.g. the inverse square potential and polynomial potentials. The numerical 
studies taking into account only one trajectory (of lowest action) show that the parameters 
of the quantum action vary smoothly with transition time T, interpolating between the 
classical action at T = and some "deep" quantum action in the asymptotic regime 
T — > oo. In the limit of large imaginary transition time, the quantum action has been 
proven to exist, and to give an exact parametrisation of transition amplitudes [18[- Also 



in this limit analytical relations exist between the the classical potential, the quantum 



potential and wave functions |l9 1 . 

What is the physical interpretation of the quantum action and its trajectories? In the 
path integral representation of the transition amplitude occurs the classical action and all 
kinds of paths (zig-zag curves). In the quantum action representation occurs the quantum 
action but only some trajectories (differentiable curves). Thus we interpret the quantum 
action as a renormalisation effect on the classical action. Quantum fluctuations appear as 
tuned parameters in the new action. This situation is similar to renormalisation in solid 
state physics: A charged particle propagating in the solid interacts with atoms. This results 
in an effective mass and charge, different from that of free propagation. 

The local quantum action and its trajectories allows to construct a quantum analogue 
phase space. From that one can compute quantum Poincare sections, Lyapunov exponents 
or the Kolmogorov- Sinai entropy like in classical nonlinear dynamics. This gives a new 
vista on quantum chaos. First results of a numerical study have been reported in Ref. 
17 1 . The quantum action has been found useful also to characterize quantum instantons 
Iq . The definition of the quantum action looks conspicously similar to that of Bohm's 
parametrisation of the wave function. So what is the difference and what is new here com- 
pared to previous studies of quantum chaos based on Bohm's approach? First, in Bohm's 
approch considers the quantum system in a state given by a particular wave function. 
Hence the fluid dynamical motion depends on such wave function (via Bohm's quantum 
potential). The quantum action, however, is the same for all initial and final states of the 
transition amplitude, and depends only on transition time. Of course, the trajectories of 
the quantum action depend on initial and final boundary conditions. Second, in Bohm's 
approach, the mass of the particle occuring in the fluid dynamics picture is the same mass 
as in the original Hamiltonian. In the quantum action, the mass is renormalized and in 
general differs in value from the classical mass. Finally, in contrast to Bohm's approach, 
the quantum action is very close in functional form to the classical action, and hence allows 
a direct and detailed comparison of phase space, as will be shown below. 



II. EXISTENCE AMD ANALYTICAL FORM OF QUANTUM ACTION IN THE 
LIMIT OF LARGE IMAGINARY TIME 



The proof of existence of the quantum action in the hmit of large imaginary time for 
Hamiltonian systems has been given for 1-D in R,ef.jl8l|. Because we consider in this work 
a 2-D Hamiltonian system, and for sake of self consistency we will outline the proof here 
for arbitrary dimension (in particular D = 2,3). Moreover, we present nonlinear differen- 
tial equations relating the classical potential to the quantum potential. Also the relation 
between the quantum action and supersymmetric quantum mechanics will be discussed. 
Let us go over to imaginary time t -^ —it. Then the transition amplitude becomes the 
Euclidean transition amplitude 



GE{xfi,T;Xin,0) = {xfi\e' 



~HT/h\ 



Xin) 



[dx] exp 



--Se[x] 



Xfr,T 



the classical action becomes the Euclidean classical action 



Se[x\ 







dt < -mx + V{x) 



and the quantum action becomes the Euclidean quantum action 



Se[x\ 



dt 



mx + V{x) 



(2) 



(3) 



(4) 
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Our goal is to prove for T — ^ oo that the Euclidean transition amplitudes can be expressed 
in terms of the Euclidean quantum action, using a single trajectory, 



Gsixfi, T; Xin, 0) = Ze exp 



Xf^,T 



•^in 5*J 



>^E[Xtraj _ 



Xf,,T 



•^in )U 



1 ~ 



dt 
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^^XffQ^j + V yXtraj J 
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(5) 



A. Necessary and sufficient conditions for existence of the quantum action 



For simplicity of notation we drop the subscript Euclidean. Let us make some assump- 
tions on the potential V{x): Let V{x) > 0. Let V{x) be a smooth (sufficiently different iable) 



function of x and let V{x) -^ oo when \x\ — > oo. Moreover, let V{x) have a single min- 
imum. Also we assume that the ground state is non-degenerate. Because the Euclidean 
transition amplitude is given by a (Wiener) path integral with a positive Euclidean action, 
G{y, T; X, 0) > for all x, y. Hence we can introduce a real function rj such that 

G{y, T; X, 0) = Go exp[-ri{y, x)] , (6) 

where Go is some constant (for fixed T) which takes care of the fact that G has a dimension 
(1/L^). 

Comparing the parametrisation of G in terms of the function rj, Eq®, with its 
parametrisation in terms of the quantum action, Eq.(0), this looks similar. In order to 
prove it we need to establish that ri{y,x) can be expressed in terms of a local action S 
evaluated along its trajectory Xtraj, such that 



V{y,x) = -S[Xtraj] 



y,t=T 

and Go = ^ . (7) 

x,t=0 



This is a necessary and sufficient condition for the existence of the quantum action. In 
order to establish that those equations hold, we proceed in the following by establishing 
a number of equations equivalent to Eq.((7)). At the end, we show that the last of those 
conditions can be satisfied. 

First, we compute partial derivatives of rj and E. Consider the functional derivative 
up to first order of the action S around the stationary point, admitting also variation of 
boundary pounts, denoted by 

~ ~ — * 

X{t) = Xtraj (t) + h{t) , 

Xtraj {t = 0) = a, Xtrajit = T) = b , 

h{t = 0) = 6a, hit = T) = 6b . (8) 



We obtain 



SS[^] = f)traj{T) ■ 5b - f)traj{0) " 5d + 0{h^) . (9) 



On the other hand, one has 



Srj{b, a) = V^ r]{y, a)\^^^ -Sb + Vg r]{b, x)y=s ■ Sa . (10) 

Comparing Eqs. ()9ll0|) for terms hnear in 6a and 6b, respectively, we find the following 
necessary and sufficient conditions for the existence of the quantum action, 

-hPtrajiO) = Vg iS|| = V^ ri{y, x) \ 

> for all boundary points {y,x) . (11) 

kPtraj{T) = Vg \t\l = Vj7 r]{y, x) J 
Eq. ffTT|) is equivalent to 

— S|| = ri{y, x) modulo a global constant . (12) 

The global constant can be absorbed into the constants Go and Z, respectively, and hence 
Eq. (fTT|) is equivalent to Eq.((7j). 



B. Use of energy conservation 

It remains to be shown how that condition Eq. ljllll can be satisfied. We do this by 
employing the principle of conservation of energy. The action given by Eq.(^ describes a 
conservative system, i.e., the force is derived from a potential and energy is conserved. In 
imaginary time, energy conservation reads 

-fkin + V = e = const , (13) 

which implies 

m - Via) = ^(vL^f - ^{f^a.r ■ (14) 

Thus combining Eq.()ll|) and Eq. ()14j) . we find another necessary and sufficient condition 
for the existance of the quantum action: The quantum action exists and is local, if there 
is a mass m and a local potential V{x), such that 



Vib) - V(a)] = (V^- rjiy, a) | .^g) ' - (v,^ r/(6, x) \s=a) ' holds for all a, b. (15) 



C. Feynman-Kac limit 

Finally, we want to show that Eq.()15|) can be satisfied in the limit T -^ oo. In this limit 
holds the Feynman-Kac formula 

G{y, T; X, 0) -^t^oo (y|V'3,.)e-^-'^/^V'p,.|x) , (16) 

where ipgr is the ground state wave function and Egr the ground state energy. Here we use 
the assumption that the ground state is not degenerate. Eq.® implies 

Goe~^^y^^^ -^T^oo (y|^3.)e-^-^/^^3,|f) . (17) 

From this we compute 

— * — * 

Vj; Viy,x)\^=b,^=d -^T-.oo= — --^ . 

4jgr{b) 



V. viy, x),.M-. — ^-- -^p^j^ ■ (18) 

Then the general condition, Eq.()15p becomes 

^|V(6) - VM] _._ (^-^-^ j - (^^-^ j for all a. i. (19) 

This holds if we can find fh and V{x), such that 
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'^^nx)-yo)= W?l for all X, (20) 



/l2 V J \ ijjgrix) 

where Vq denotes the minimum of the potential. This condition can be satisfied. This 
establishes the existence of a local quantum action and finishes the proof. 

D. Relation between classical and quantum potential 

We continue to work in the limit T — ^ oo. Let us consider the stationary Schrodinger 
equation for the ground state, which can be written in the form 

Aipgrlx) 2m 



tpgr{x) h? 



\V{x)-Eg,], (21) 



where m is the classical mass and V is the classical potential. Because the transition 
amplitude is positive, also the ground state wave function obeys ipgr{x) > for all x. 
Hence we can define the function U{x) 

Uix) = \og^Pgrix) . (22) 



Using Eqs. ()20|21|) we compute 



AU{x) = ^[\/(f) - Vo] + |^[\/(f) - E,r] . (23) 



It can be cast into the form 

2171 

■) + {yu{x)y = 

2m ,~ 



AU{x) + {VU{x)y = -^[V{x) - Egr] , 



^[Vix)-Vo] = i\/Ur. (24) 

The first equation is a generalized Riccati-type differential equation for the function U{x). 
From its solution one can obtain the the quantum potential, more precisely mV{x). Hence, 
Eqs. (l24|) constitute a transformation rule 

m,V{x) — > ThV{x) . (25) 



From the point of view of numerical solution, it may be advantageous to cast Eq. ()24j) in a 
slightly different form. Let us define 

W{x) = VU{x) . (26) 

Then the generalized Riccati differential equation becomes 



^gr] 



— * — * — * ^TD 

V ■ W{x) + {W{x)Y = -^[V{x) - E,^ 

Off} ^ ^ —^ 

— [V{x)-V,] = {W{x)f. (27) 



E. Relation between quantum action and potentials in supersymmetric quantum 
mechanics 



starting with Witten. wo^^ it wa= recognized that snpe.y„nnet.y can be appUed to 

quantum mechanics. In Ref.|22| the correspondence between the 1-D Fokker-Planck equa- 
tion with a bi-stable potential and supersymmetric quantum mechanics has been shown to 
facihtate the calculation of a small eigenvalue that controls the rate at which equilibrium is 
approached. Kumar et al. ^233 used the supersymmetric approach to compute the tunneling 
exchange integral t = Ei — Eq and pointed out the Riccati differential equation between 
the classical and supersymmetric partner potential. Let us consider a 1-D Hamiltonian 
system, which in the notation of Ref . J2J] and using h = 2m = 1 is given by 

H^ = -^,+V.{x), (28) 

where V- denotes the standard classical potential. One can construct a supersymmetric 
partner potential V^{x), 

V^ix) = V4x) - 2-f ^ = -V4x) + 2(^)^ . (29) 

dX Ipgr 'Ipgr 



Defining a superpotential Ws{x) by 



one obtains 



Ws{x) = -^ , (30) 

Ygr 



V^{x) = W^ix) ± ^ , (31) 

which is a Riccati-type equation. The following similarities are apparent. Comparing 
Eqs.dlSEnD one finds 

Wix) = -Wsix) , (32) 

which satisfy the corresponding Riccati differential equations (J27l3ip . The following re- 
marks are in order. In the super-symmetric approach the mass is identical to the classical 
mass, while in the quantum action, the mass gets tuned, in general. In contrast to the 
super-symmetric quantum Hamiltonian, the quantum action gives a parametrisation of 
quantum transition amplitudes (see Eq.((T))), which requires a tuning of the mass. 
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III. MODEL AND ITS QUANTUM ACTION 
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FIG. 1: Ground state wave function iljgr{x^y). 



In order to study full chaos the so-called K-system (2-D Hamiltonian, potential V — 
x^y'^) is widely used. It is almost globally chaotic, having small islands of stability J25|. 
In order to study mixed dynamics (entangling chaos and regular islands) it is numerically 
convenient to consider the following classical system |26i |. 

r-T ^ 

S= I dt -m{x^ + y") - V{x, y) , V = V2{x^ + y^) + t)22xV • (33) 
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FIG. 2: Quantum potential rh{V{x,y) — Vmin)- 

The parameters of the classical action are given in Tab.p. The ground state wave function 
of the quantum system, obtained by solving the Schrodinger equation, is depicted in Fig. [l]. 
The quantum action (see Eq.((H)) takes the form 



S 



dt -m{x^ + ij^) — V{x, y) 



(34) 



We need to determine the quantum mass (tuned mass) fh and the quantum potential 
V{x,y). For large imaginary time T we have analytical results on the form of V. The 
function m(V^(x, y) — Vmin) is given by Eq. fj^Uj) or equivalently by the solution of Eq.(j2Zj). 
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TABLE I: Parameters of classical action vs. quantum action. ^22 = 0.25, T = 4.5. 



Parameter 


Class. Action 


Quant. Action 


m 


1 


0.976(8) 


Vo 





1.3992(4) 


V2 


0.5 


0.5684(3) 


■"22 


0.25 


0.2469(4) 


va 





-0.00067(7) 



This is depicted in Fig. j2j . We need to determine m and V separately and also the parameter 
Z of the quantum action. Those have been computed numerically in a non-perturbative 
way from a global fit to transition amplitudes in imaginary time. For this purpose we made 
a polynomial ansatz for the quantum action of the following form, 



V{x,y) = vq + viixy + V2{x + y ) + V22X y + vi^ixy + x y) + V4{x +y) 



+ ^^24(3; y +x y ) +vux y 



(35) 



We have computed transition amplitudes G{xfi,tfi;Xin,tin) by solving the Schrodinger 
equation in imaginary time keeping the transition time T fixed. The determination of un- 
known parameters of the quantum action, like Z, rh, etc., requires as many equations as 
unknowns. In practice we use much more equations (overdetermined system), by consider- 
ing a large set of pairs of initial and final boundary points {xj,Xi) and the corresponding 
transition amplitudes. Next we make an initial guess for the unknown parameters, say Z^'^^ 
ifi^^\ etc. and solve the Euler-Lagrange equations of motion derived from the quantum ac- 
tion S^^\ for each pair of boundary points. Then we compute the error 



N 



G{xi, T; £,-, 0) — Z^^' exp 



1 c(o)rHo) 1 



Xi,T' 



x,,0 



(36) 



In order to find eventually the physical values of the unknown parameters we apply a 
variational principle, i.e. we search in the multi-parameter space to minimize the error 
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globally (over all pairs of boundary points). For more details see Ref.|l5i|. As a results, 




2000 4000 6000 8000 10000 12000 14000 16000 18000 20000 
Temps 



FIG. 3: Typical time evolution of Lyapunov exponents. ^22 = 0.25. E = 2. 

the terms Vu, V13, V4, V24, v^i, were found to be quite small or compatible with zero. 
The parameters of the corresponding quantum action for transition time T = 4.5 (large 
compared to dynamical time scale Tsc = l/Egj.) and parameter f22 = 0.25 (which controls 
chaos) are shown in Tab.p. The bracket gives the estimated error of the fit. The numerical 
results presented below refer to ^22 = 0.25 and ^22 = 0.05. Before we present our 
comparison of classical to quantum chaos, we ask what behavior do we expect? In the 1-D 
quartic potential quantum effects produce a strong positive quadratic term in the quantum 
action [1^. In the 1-D double well potential, the quantum potential becomes softened, 
i.e., it has a lower barrier and closer minima. As a consequence, the instanton solutions of 
the quantum action are softer than the corresponding classical instantons [16]. Are such 
softening effects also occuring in chaotic phenomena? 

IV. CLASSICAL CHAOS VS. QUANTUM CHAOS 



The notion of chaotic versus regular phase space is based on Lyapunov exponents (A > 
versus A < 0). Numerically it is difficult to distinguish a small positive from a zero 
Lyapunov exponent. In order to discriminate those cases we followed the time evolution for 
a long time. We used Tc = 20000 to measure A. This is depicted in Fig. |Hj. A comparison of 
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FIG. 4: Poincare section 
E = 2. 



of classical system (top) and quantum system (bottom). ^22 = 0.25. 
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FIG. 5: Poincare section of classical system (left) and quantum system (right). ^22 = 0.25. E = 6. 

Poincare sections of the classical and the quantum system is shown in Fig. ji] for ^22 = 0.25 
and E = 2. One observes a smaller number of chaotic trajectories in the quantum case than 
in the classical case. In particular, this is visible near to the hyperbolic fixed points (see 
insert). The same comparison for some higher energy {E = 6) is shown in Fig.jS]. With 
increasing energy the difference between classical and quantum phase space becomes more 
pronounced. In order to get a quantitative measure for chaotic versus regular phase space, 
we have chosen randomly a large number of initial conditions and computed the finite time 
Lyapunov exponent for each such trajectory. The distribution of these Lyapunov exponents 
in the neighborhood of A = is shown in Fig. [0], and those for < A < 0.25 is shown in 
Fig.[7|. Fig.jU] shows the distribution of Lyapunov exponents for different energies for the 
classical and quantum system. While A = corresponds to regular behavior, A > indicates 
chaotic behavior. For the purpose to distinguish numerically the two regimes, we used this 
distribution to define some cut-off Ac (vertical dashed line in Fig. [Hj). Fig.jBj also shows the 
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FIG. 6: Distribution of Lyapunov exponents near A = for different energies E. Classical system 
(left column) and quantum system (right column). W22 = 0.25. 

cumulative distribution P(A) (full line), which is a measure for the degree of regularity. As 
can be seen this curve is higher for the quantum systen compared to the classical system. 
Moreover, Fig. |B] displays the coefficient R which denotes the ratio of chaotic phase space 
over total phase space. The distribution of positive Lyapunov exponents is displayed in 
Fig. [21. First, one observes a pronounced peak near zero, which has been magnified in 
Fig. p. Second, with increase of energy, the system becomes more chaotic, i.e. R increases. 
Third, with increasing energy, the distribution of positive Lyapunov exponents seems to 
approach a Gaussian shape. 

Fourth, the width of the Gaussian (variance) diminishes with increasing energy. Fifth, 
with increasing energy, a strong linear correlation develops between the expectation value of 
A and energy E, which can be represented by a linear fit < A >= A*-"-' + e-E'. (fit parameters 
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FIG. 7: Distribution of positive Lyapunov exponents. Classical system (left column) and quantum 
system (right column). z;22 = 0.25. 



,(0) 



,(0) 



are: A^^' = -0.040, e^ = 0.033 and A^';^ = -0.026, eg„ = 0.024, i.e. e,^ < e^). All 



those features hold for both, the classical system and the quantum system. A quantitative 
difference between classical and quantum system is apparent in the degree of chaoticity R, 
shown in Figs. |6l7] . The corresponding results for ^22 = 0.05 are given in Figs. |8l9] . which 
show qualitatively the same behavior. The basic difference between the data for ^22 = 0.25 
and f22 = 0.05 shows up in the degree of chaoticity R which is smaller for the latter, and 
also in the location of the peak in the distribution of positive Lyapunov expenents (compare 
Figs. [30]), which is also smaller for the latter. This shows that the parameter ^22 which 
controls the degree of chaoticity in the classical system, has a counterpart in the parameter 
■y22 in the quantum system, which plays a similar role. For all energies, and all values of f 22, 
R plotted in Fig.^U] is smaller in the quantum system than in the corresponding classical 
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FIG. 8: Distribution of Lyapunov exponents near A = for different energies E. Classical system 
(left column) and quantum system (right column). W22 = 0.05. 

system, i.e. the classical system turns out to be more chaotic than the quantum system. 
This is the main result of this work. 

How can we understand the observed behavior? The distribution of Lyapunov exponents 
approaching a Gaussian with increasing energy is well known from Hamiltonian systems 
with mixed dynamics psf]. The shape of the curve and its decreasing width can be under- 
stood as a consequence of the central limit theorem, which predicts such behavior for an 
arithmetic mean of A^ random variables for N —>■ 00. On the other hand, regular islands in 
an otherwise chaotic phase space are known to give non-Gaussian tails. What can we say 
about those tails? In Hamiltonian chaotic systems, when increasing the energy then phase 
space becomes more chaotic, i.e. regular islands become smaller or disappear. Thus it is 
plausible that the non-Gaussian tails of the distribution of positive Lyaponov exponents. 
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FIG. 9: Distribution of positive Lyapunov exponents. Classical system (left column) and quantum 
system (right column). W22 = 0.05. 

associated with the regular islands, diminish. We have estimated numerically the weight of 
non-Gaussian tails using a fit of the data by a Gaussian distribution. For example, in the 
distribution of classical Lyapunov exponents corresponding to ^22 = 0.25 (see left column 
of Fig.jTj), the relative error e of such fits is as follows: E = 2, e = 0.60, E = 4, e = 0.39, 
E = 6, e = 0.30, E = 8, e = 0.27, which diminishes with increasing energy. Qualitatively 
similar behavior is found in the corresponding quantum Lyapunov exponents as well as for 
the data corresponding the coupling parameter ^22 = 0.05. 

Second, quantum fluctuations which occur in the parameters of the quantum action can 
be viewed as renormalisation effects of the classical parameters [19| . The value of f 2 (if sole 
term of potential would render the system integrable) is increased, while f 22 (term which 
drives the system away from integrability) is decreased (see Tab.jl]. Hence chaos is reduced 
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FIG. 10: Volume of chaotic phase space (positive Lyapunov exponent) over total phase space 
(hnes are guides to the eye). Left ^22 = 0.05, right 1^22 = 0.25. 

in the quantum system. Such increase of the quadratic term has been observed in other 
systems also and seems to occur more generally. 

V. CONCLUSION 

What is the significance of our findings? In particular, what is the physical significance 
of chaotic properties of trajectories for the quantum system? (i) We found that the quan- 
tum action for large transition times allows to construct a quantum analogue phase space 
for Hamiltonian systems with mixed chaotic dynamics. Such phase space resembles the 
classical phase space, in particular, when the parameter driving chaos is small. For the 
mixed system considered, the quantum analogue phase space, like classical phase space, 
shows islands of regular behavior interwoven with chaotic regions. Moreover, hyperbolic 
fixed points occuring in classical phase space survive in the quantum phase space. We 
compared quantitatively the classical system with the quantum system and found that the 
latter is systematically less chaotic than the former. The reason for such behavior, in our 
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opinion, is the tendency of quantum fluctuations to drive the quantum action nearer to a 
Gaussian fixed point. Moreover, we have investigated the quantum analogue phase space as 
a function of energy. We find that such phase space, like the classical phase space, becomes 
more chaotic when energy increases. Thus introducing a quantum analogue phase space 
gives a detailed and quantitative measure reflecting how much chaos grows as a function 
of energy. 

(ii) The significance of trajectories can be seen by drawing an analogy with the in- 
terpretation of scattering experiments. In a scattering experiment, one usually measures 
counting rates in detectors, which when positioned at different angles give a differential 
cross section. Using theoretical input, like the S-matrix, unitarity, and symmetries like 
rotational symmetry of the potential, one can obtain derived quantities, like scattering 
phase shifts Si{E). They give a picture of scattering as function of energy and indicate 
the existence of resonances in certain channels (jumps of phase shift). The partial wave 
expansion in conjunction with phase shifts Si{E) gives a representation of the S-matrix and 
of scattering cross sections. In analogy to this picture, we consider the quantum action 
and its trajectories as a parametrisation of the quantum mechanical transition amplitudes. 
The quantum action and its trajectories are derived quantities. These trajectories in con- 
junction with tools of nonlinear dynamics (Poincare sections, Lyapunov exponents) give a 
picture of (makes "visible") the chaotic dynamics in the quantum system. One may view 
the Lyapunov exponents computed from trajectories of the quantum action in analogy to 
scattering phase shifts computed from a partial wave analysis of scattering cross sections. 

(iii) The trajectories of the quantum analogue phase space gain a particular importance 
in the interpretation of of the phenomenon of chaos assisted quantum dynamical tunneling. 
This phenomenon is similar to tunneling through a potential barrier, however, in this case 
the classical transition is forbidden by a dynamical law, the presence of a KAM surface. 
The experiments with ultra cold atoms in an amplitude modulated optical standing wave 
reported by Hensinger et al. |2^ and Steck et al. (2i| demonstrate that such dynamical 
tunneling exists and its amplitude is enhanced by chaos. The tunneling occurs between 
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regions which in classical phase space correspond to regions of stable motion (fixed points). 
In the quantum system this region is associated with eigenstates of the Floquet operator. 
Hensinger et al. point out that "Floquet states do not necessarily overlap well with the 
regions of regular motion of the classical Poincare map, except in the semi-classical limit" . 
Steck et al. work with wave packets narrow in momentum space but not localized in po- 
sition. In order to get a picture of phase space they use a classical phase space with a 
distribution with the same position and momentum marginal distribution as the Wigner 
function. We suggest that in order to better understand the mechanism of dynamical tun- 
neling, the role of localized states, decoherence, and a higher number of Floquet states, a 
quantum analogue phase space and its trajectories from the quantum action are a com- 
plementary tool and should give additional useful insight in the dynamics of the quantum 
system. 
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